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TECHNICAL MEMORANDUM 


PRELIMINARY CONTROL SYSTEM DESIGN AND ANALYSIS FOR THE 
SPACE STATION FURNACE FACILITY THERMAL 
CONTROL SYSTEM 


INTRODUCTION 


The Space Station Furnace Facility (SSFF) is a facility for materials research in the microgravity 
environment of the international Space Station Alpha United States Laboratory. The SSFF provides the 
necessary core systems to operate various materials processing furnaces. The thermal control system 
(TCS) is defined as one of the core systems, and its function is to collect excess heat from furnaces and 
to provide precise cold temperature control of components and of certain furnace zones. The TCS 
utilizes single-phase water as its cooling medium and consists of piping, heat exchangers, cold plates, 
valves, flow sensors, temperature sensors, pressure sensors, and a pump. The TCS design layout is 
shown in figure 1. Combinations of the subsystems are plumbed in parallel to help heat rejection con- 
trollability. 

The specific control system objective for the TCS is to ( 1 ) control the total pressure drop across 
the parallel system and (2) track defined flow rate profiles in each parallel subsystem path. Two control 
inputs are employed to achieve the control system objective. One control input is the pump speed com- 
mand, which is used to modulate pump speed to maintain the desired system pressure drop. This serves 
the purpose of decoupling the flow dynamics between the respective parallel paths. The other control 
input is the valve speed command, which is used to modulate the flow control valve to track defined 
flow rate profiles in each subsystem path. Tracking the desired flow profiles assures that subsystem tem- 
peratures are maintained. 


PURPOSE AND SCOPE OF ANALYSIS 


The SSFF TCS control system analysis was performed to justify and fine-tune digital control 
algorithms that generate the control inputs. The analysis is an ongoing activity consisting of model 
development, model validation, and open- and closed-loop analyses utilizing frequency, time, and dis- 
crete domain techniques. This report serves as the documentation of the TCS control system design and 
analysis performed prior to the SSFF preliminary design review. The mathematical modeling, control 
system designs, and control system analyses were performed using Marshall System for Aerospace 
Simulation (MARSYAS), a control system analysis and time simulation package that was developed by 
the Control Systems Division at Marshall Space Flight Center. The following sections are included in 
this report: the description of the TCS mathematical model; the control system design, both continuous 
and discrete; the control system analysis including stability margins; the system performance analysis 
including time simulations; and the MARSYAS code used in the simulations and analyses. 



CONTROL SYSTEM MATHEMATICAL MODEL 


The TCS has a pump source that forces cold water through six subsystem flow paths (fig. 1 
shows eight paths, but only six are operated simultaneously) plumbed in parallel. The water is used to 
collect waste energy and to control thermal environments of the various subsystems. The mathematical 
model of the TCS is divided into two submodels: the flow path model and the source model. 


Flow Path Model 

Each flow path consists of components (cold plates, heat exchangers, flex lines, quick discon- 
nects, furnace cooling jackets, etc.) that are plumbed in series with a flow control actuator and valve. 
The flow model has an equation that describes the mass flow rate of water and two equations that 
describe the flow control actuator and valve. The mass flow rate equation (1) is derived from an equiva- 
lent circuit (fig. 2) analogy using lumped parameters for flow resistance, inertance, and capacitance. 1 
The flow equation is written as: 


dW L _j_ 
dt ~ L , 
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where W, is the mass flow rate through each path, 0, is the angle of the flow control valve in each path, 
6% is the velocity of the flow control valve in each path, L,- is the lumped flow inertance through each 
path, a , ft, and 5 are the valve resistance constants, /?, is the lumped flow resistance constant for each 
path, and P is the system pressure drop across the paths. 

The flow control actuator consists of an analog speed controller and a brushless dc motor. The 
flow control actuator and valve are connected through a shaft and together provide precise control of the 
flow profiles in each path. The valve equation (2) is derived from the idea that the valve acts as an ideal 
integrator of its commanded speed. The valve equation is written as: 


UU: w , ^ 

= - Vi = 1,6 . (2) 

The speed controller equation (3) is approximated from the idea that the controller achieves the 
commanded valve speed with a slight time delay. The speed controller equation is written as: 

do), 1 

-jf = jr(k v u r a) i ) , Vi =1,6, (3) 

where ris the time constant delay of the commanded valve speed, k v is the conversion constant between 
voltage and valve speed, and w, is the valve input voltage command for each path. 


Source Model 

The source model consists of the system flow and pressure drop, pump motor, pump speed con- 
troller, and system pressure controller equations. The system flow equation (4) and the system pressure 
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drop equation (5) are derived using an equivalent circuit analogy. 1 The system flow equation is written 
as: 


dW T i 
dt ~ L eq 


-R eq Wl+Pp(W T ,Q)p)+P , 


(4) 


where W T is the system flow through the pump, P is the system pressure drop across the paths, L eq is the 
equivalent system flow inertance, R eq is the equivalent system flow resistance, ci)p is the speed of the 
pump, and P p (W T ,C0p ) is the pump map function (fig. 3). The system pressure drop equation is written as: 


dP _ 2 
dt C eq 



(5) 


where C eq is the equivalent system flow capacitance, and Wj is the flow through each subsystem path. 

The dynamic model representation of the pump and the corresponding parameter data are 
obtained from the pump vendor, Allied Signal. 2 The pump motor and controller equations are written as 

follows: 
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- 3 JL = FX c+ G [o fbk P P set ] T , (9) 

where o p is the pump speed, I a is the pump motor armature current, a>p, k is the pump speed feedback, X c 
is the 2x1 row-vector of controller states ,j m i is the rotor plus load inertia, k; is the current to motor 
torque gain, A: i is the speed to load torque gain, hm is the pump motor inductance, R m is the pump motor 
resistance, k e is the motor back emf gain, k 2 is the controller feed forward gain, k 3 is the controller feed 
forward gain, E is the 1x2 row-vector of controller gains, is the speed filter time constant, F is the 2x2 
matrix of controller gains, and G is the 2x3 plant state input matrix. 


CONTROL SYSTEM DESIGN 


The specific control system design objective is to (1) control the total pressure drop across the 
parallel system and (2) track defined flow rate profiles in each parallel subsystem path. Two controllers 
are designed to achieve this objective. The first is a pressure controller that is used to modulate the pump 
speed to maintain the total pressure drop. The second is a flow controller that is used to modulate the 
flow control valve to track defined flow profiles in each subsystem path. 
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Pressure Controller 


The coupling effect between subsystem paths is the most significant disturbance to controlling 
flow. For instance, flow rate modulation in one subsystem path causes a disturbance to the other subsys- 
tem paths. This so-called coupling effect can be seen mathematically from the plant equations. When 
written in matrix notation, the plant’s flow equations have off-diagonal terms in pressure. These terms 
represent the coupling between subsystem paths. To minimize this effect, the pressure between the sub- 
system paths is controlled to a constant. The continuous and discrete control system designs to perform 
this pressure control are shown in Figure 4. The continuous system was designed and then converted to a 
discrete equivalent system 3 with a sample rate of 1 s. 

The pressure control system design has a lead-lag filter that improves system stability and 
response through the placement of poles and zeros in the frequency domain. The prime purpose of a lead 
filter is to add phase near the crossover frequency in order to increase the phase margin. Accompanied 
with this phase lead is a gain increase that will increase the crossover frequency. The basic idea behind 
the lag filter design is to reduce the gain at “middle” frequencies in order to reduce the crossover fre- 
quency to a lower value than for the uncompensated system. Combining the lead and lag filters provides 
attenuation below the crossover frequency to decrease the phase lag at crossover and phase lead closer to 
the crossover frequency in order to increase the phase lead of the uncompensated system at the crossover 
frequency. 4 

The output of the controller is the pump speed command that is limited and held for a 1-s sample 
delay. The MARSYAS code implementing the continuous and discrete controllers is listed in the 
appendix. 


Flow Controller 

Performance requirements for flow control such as rise times, overshoots, and settling times are 
not specified for a good reason — they are not critical. Because the dynamic lag between flows and tem- 
peratures is so large (minutes), transient variations in flow and pressure do not significantly affect sub- 
system temperatures. However, long-term variations in flow do affect temperatures. With this in mind, 
the control goal is to design a very stable control system that is insensitive to disturbances and has good 
steady-state accuracy. 

The continuous and discrete control system designs for flow rate control to each EM and to the 
core rack are shown in figure 5. The continuous system was designed and then converted to a discrete 
equivalent system 3 with a sample rate of 1 s. The flow set point and the flow error signals are low pass 
filtered to smooth out the system’s response to disturbances. The flow control system design has a pro- 
portional-plus-derivative (PD) control algorithm 5 that achieves the desired flow tracking performance. 
The valve model has a free integrator, which coupled with the PD controller, achieves zero steady-state 
flow tracking error. A valve command limiter and a computer processing delay are also included. The 
MARSYAS code implementing the continuous and discrete controllers is listed in the appendix. 


CONTROL SYSTEM ANALYSIS 


An analysis is performed to determine the system’s sensitivity to variations in pressure and flow. 
The control system shown in figure 4 is used to analyze the system’s sensitivity to pressure variations. 
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The feedback control loop is “opened at” the Err variable and a new input and output created. The trans- 
fer function between the new input and output is generated and a frequency response run. The response 
is shown in figure 6. From this plot, the gain and phase margins are obtained as 20 dB and 1 10 , respec- 
tively. 

The gain margin is defined as the amount of gain the system would require to cross 0 dB at 180 
of phase. The phase margin is defined as the amount of phase the system would require to cross 180 at 
0 dB of gain. Gain and phase margins quantify the relative stability or sensitivity that the system has to a 

certain parameter. 


The control system shown in figure 5 is used to analyze the system’s sensitivity to variations in 
flow. The feedback control loop is “opened at” the Err variable and a new input and output created. This 
process is repeated for all three flow control loops (two EM’s and core). The resulting frequency 
response, shown in figure 7, is representative of all three flow control systems (EMI A, EM2B and the 
core). The plot shows a gain and phase margin of 50 dB and 78°, respectively. These large stability 
margins indicate that the system is very stable to variations in flow and pressure. 


CONTROL SYSTEM SIMULATIONS 

The integrated TCS is simulated on the digital computer. Active control is delivered to the sys- 
tem at a sample rate of 1 s. For this discussion, the pressure and flow nomenclature shown in figure 2 is 
used For this simulation, the control objective is to track flow rate profiles for W2, W13, and W18, 
which represent the flows through EMI A, the core, and EM2B, respectively, and to maintain the static 
flows for W3, W6, and W17, which represent the flows through the constant flow subsystem paths. io 
help decouple the flow dynamics between the subsystem paths, the differential pressure between PI and 
P10 is controlled to 10 lb/in 2 . 

The model is initialized at a steady-state operating point and simulated for 6 min. The input to 
the model is simultaneous changes in all three flow rate profiles for W2, W13, and W18. The solid and 
dotted lines on the top graphs in figures 8, 9, and 10 are plots of the flows and profiles, respectively. The 
bottom graphs in figures 8, 9, and 10 are plots of the respective valve positions that are modulated by the 
flow control algorithm to track the flow profiles. The top graph in figure 1 1 is the plot of the controlled 
system pressure drop. The bottom graph in figure 1 1 is the commanded pump speed to maintain the 
pressure drop. The graph in figure 12 is a plot of the total system flow through the pump, and the graphs 
in figure 13 are the plots of the static flows in the constant flow subsystem paths. 


CONCLUSIONS 

The nonlinear mathematical model of the TCS has been developed and presented in this report. A 
very robust control system has been designed to maintain flow profiles and pressures in the TCS in the 
face of simultaneous disturbances. Operationally, the TCS will not experience the simultaneous high- 
frequency flow profile changes applied in the simulation; nevertheless, these disturbances were applied 
to assess its performance in a worst-case scenario. The simulation demonstrates that the design goal of a 
very stable control system that is insensitive to disturbances and has good steady-state accuracy is 
achieved. Control system analyses and simulation runs show that the system is very stable and has satis- 
factory transient and steady-state response. The control system designs and analyses will continue to 
mature as model validations are performed. 
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Continuous Pressure Controller 



Discrete Pressure Controller 



Figure 4. Pressure controller block diagrams. 


Continuous Flow Controller 



Discrete Flow Controller 



Figure 5. Flow controller block diagrams. 
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Figure 6. Pressure frequency response. 




Figure 7. Flow frequency response. 
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Figure 8. EMI A flow and valve angle plots. 
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Figure 9. Core flow and valve angle plots. 


10 







w 

1 

8 

L 

B 

m 

/ 

H 

R 




TIME (MINS) 

Figure 10. EM2B flow and valve angle plots. 
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Figure 11. Pressure drop and pump speed plots. 
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APPENDIX 


**** main **** 

* This file is located in /home/jacksme/marsyas/furn/main on host edl4j. 

* It contains the main program for the SSFF TCS simulation. It also 

* contains defined constants and the respective inputs and outputs of 

* the main model and the various submodels $ 

* program global constant definitions $ 

CONSTANT: PI = 3.1415927 $ 

* density of water (lbm/ft A 3) $ 

: RHOF = 62.268 $ 

* conversion constant (lbm-ft / sec A 2-lbf) $ 

:GC = 32.174$ 

* inner diameter for components (ft) $ 

: DC = 0.305/12 $ 

* inner diameter for pipes (ft) $ 

:DI = 0.5/12$ 

* algebraic loop breaker delay time (secs) $ 

: TDELAY = 0.005 $ 

* integration timestep (secs) $ 

: TSTEP = 0.001 $ 

* simulation stoptime (secs) $ 

: STOPT = 0.0 $ 

* include the file /home/jacksme/marsyas/furn/defines which has all 

* function call definitions $ 

INCLUDE: "defines" $ 

* main model definition $ 

* inputs: 

* amfl - valve angle for static path in irl (degs) $ 

* amc - valve angle for static path in core (degs) $ 

* amf2 - valve angle for static path in ir2 (degs) $ 

* pump\ref\speed - 100 percent pump speed (rpm) $ 

* t - water temperature (degs f) $ 

* pheadset - system pressure set point (psi) $ 

* w2set - flow set point through em la in irl (lbm/hr) $ 

* wl3set - flow set point through core (lbm/hr) $ 

* wl8set - flow set point through em2b in ir2 (lbm/hr) $ 

* output definitions: 

* time\min - simulation time (mins) $ 

* phead - system pressure (psi) $ 

* pump\speed - pump speed (rpm) $ 

* w2set - flow set point through emla in irl (lbm/hr) $ 

* w2 - flow through emla in irl (lbm/hr) $ 

* emla\valve\angle - valve position in path emla of irl (degs) $ 

* wl3set - flow set point through core (lbm/hr) $ 

* wl3 - flow through core (lbm/hr) $ 
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* core\valve\angle - valve position in core (degs) $ 

* wl8set - flow set point through em2b in ir2 (lbm/hr) $ 

* wl 8 - flow through em2b in ir2 (lbm/hr) $ 

* em2b\valve\angle - valve position in path em2b of ir2 (degs) $ 

* wl - total system flow through pump (lbm/hr) $ 

* w3 - flow through static path in irl (lbm/hr) $ 

* w6 - flow through static path in core (lbm/hr) $ 

* wl7 - flow through static path in ir2 (lbm/hr) $ 

MODEL: MAIN $ 

INPUTS: AMFl,AMC,AMF2,PUMP\REF\SPEED,f,PHEADSET,W2SET,W13SET W18SET $ 
OUTPUTS: TIME\MIN,PHEAD,PUMP\SPEED,W2SET,W2,EM 1 AWALVEVANGLE.W 1 3SET W 1 
CORE\VALVE\ANGLE,W 1 8SET.W 1 8,EM2B\VALVE\ANGLE,W 1 , W3.W6.W 17$ 

* pump speed controller submodel definition $ 

* inputs: 

* pump\speed\cmd - pump speed command from pressure controller (rads/sec) $ 

* pump\speed\sen - pump speed (rads/sec) $ 

* outputs: 

* pump\volt\cmd - pump motor voltage command (volts) $ 

SUBMODEL: PUMP\SPEED\CONTROLLER; 

INPUTS: PUMP\SPEED\CMD,PUMP\SPEED\SEN- 
OUTPUTS: PUMPWOLT\CMD $ 

* pump motor submodel definition $ 

* input: 

* pump\volt\cmd - pump motor voltage command (volts) $ 

* outputs: 

* pump\speed - pump speed (rpm) $ 

* pump\speed\sen - pump speed (rads/sec) $ 

* pump\cun\sen - pump motor armature currenct (amps) $ 

SUBMODEL: PUMP\MOTOR; 

INPUTS: PUMPWOURCMD; 

OUTPUTS: PUMP\SPEED,PUMP\SPEED\SEN,PUMP\CURR\SEN $ 

* flow control valve\actuator submodel definitions $ 

* inputs: 

* 'patlAvalveXcmd - valve command (volts) $ 

* outputs: 

* ’patlAvalveWigle - valve angle (degs) $ 

SUBMODEL: EMI A\VALVE\ACTUATOR ; 

INPUTS: EM 1 AW ALVE\CMD ; 

OUTPUTS: EM 1 AW ALVE\ANGLE $ 

SUBMODEL: CORE\ V AL V E\ ACTU ATOR ; 

INPUTS: CORE\VALVE\CMD ; 

OUTPUTS: CORE\VALVE\ANGLE $ 

SUBMODEL: EM2B\VALVE\ACTUATOR ; 

INPUTS: EM2B\VALVE\CMD ; 

OUTPUTS: EM2B\VALVE\ANGLE $ 
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* continuous and discrete pressure controller submodel definitions $ 

* inputs: 

* phead - system pressure (psi) $ 

* pheadset - system pressure setpoint (psi) $ 

* outputs: 

* pump\speed\cmd - pump speed command (rads/sec) $ 
SUBMODEL: PRESSURE\CONTROLLER\CONTINUOUS; 
INPUTS: PHEAD.PHEADSET ; 

OUTPUTS: PUMP\SPEED\CMD $ 

SUBMODEL: PRESS URE\CONTROLLER\DI SCRETE; 

INPUTS: PHE AD.PHE ADSET ; 

OUTPUTS: PUMP\SPEED\CMD $ 

* continuous and discrete flow controller submodel definitions$ 

* inputs: 

* w# - path flow (lbm/hr) $ 

* w#set - path flow set point (lbm/hr) $ 

* outputs: 

* 'patlAvalve\cmd - path valve command (volts) $ 

SUBMODEL: EMlA\FLOW\CONTROLLER\CONTINUOUS ; 
INPUTS: W2.W2SET ; 

OUTPUTS: EM1A\VALVE\CMD $ 

SUBMODEL: EMlA\FLOW\CONTROLLER\DISCRETE ; 
INPUTS: W2.W2SET ; 

OUTPUTS: EM1A\VALVE\CMD $ 

SUBMODEL: CORE\FLOW\CONTROLLER\CONTINUOUS ; 
INPUTS: W 13, W13SET ; 

OUTPUTS: COREYV ALVE\CMD $ 

SUBMODEL: CORE\FLOW\CONTROLLER\DISCRETE ; 
INPUTS: W 1 3, W13SET; 

OUTPUTS: COREWALVEVCMD $ 

SUBMODEL: EM2B\FLOW\CONTROLLER\CONTINUOUS ; 
INPUTS: W18/W18SET ; 

OUTPUTS: EM2B\VALVE\CMD $ 

SUBMODEL: EM2B\FLOW\CONTROLLER\DISCRETE ; 
INPUTS: W18,W18SET ; 

OUTPUTS: EM2B\VALVE\CMD $ 


**** flow **** 

* This file is located in /home/jacksme/marsyas/fum/flow on host edl4j. It 

* is part of the main model and contains the flow and pressure dynamic 

* equations used in the SSFF TCS simulation $ 
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* The following equations are written utilizing kirchoffs voltage law 

* and an equivalent electrical circuit analogy where pressure and flow are 

* analogous to voltage and current, respectively. The circuit drawing used 

* to obtain these equations is shown in the TECS Data Book (MSFC-DOC-2383) 

* in the controls analysis section. Certain paths (that is, w2, w3, w6, wl3, 

* wl7, and wl8) have been modeled without flow inertance because a very large 

* eigenvalue results from the path's R/L combination. This causes the 

* simulation to be run at a very slow timestep. The flow inertances associated 

* with those paths are included in the upstream or downstream paths, hence 

* preserving the system inertance. An algebraic loop results in the no- 

* inertance paths between the calculation of flow and its resistance. 

* Algebraic loop breakers are added below to decrease simulation time. $ 

* flow equations (lbm/hr) $ 

EQUATION: Wl' = 1/L1*(PHP - PI - PL1) $ 

: W2D = (P3 - P4)/(RVF 1 A+R2) $ 

: W2' = (W2D - W2)/(TDELAY) $ 

: W3D = (P3 - P4)/(RMF1+R3) $ 

: W3’ = (W3D - W3)/(TDELAY) $ 

: W4’ = 1/L4*(P2 - P3 - PL4) $ 

: W5' = 1/L5*(P4 - P5 - PL5) $ 

: W6D = (P2 - P5)/(RMC+R6) $ 

: W6' = (W6D - W6)/(TDELAY) $ 

: W7' = 1/L7*(P1 - P2 - PL7) $ 

: W8‘ = 1/L8*(P1 - P6 - PL8) $ 

: W9’ = 1/L9*(P5 - P9 - PL9) $ 

: W10’ = 1/L10*(P6 - P7 - PL10) $ 

: Wll' = 1/Ll 1*(P7 - P8 - PL1 1) $ 

: W12’ = 1/L12*(P7 - P8 - PL12) $ 

: W13D = (P8 - P9)/(R13+RVC) $ 

: W13' = (W13D - W13)/(TDELAY) $ 

: W14’ = 1/L14*(P9 - P10 - PL14) $ 

: W15' = 1/L15*(P6 - Pll -PL15)$ 

: W16’ = 1/L16*(P12 - P10 - PL16) $ 

: W17D = (Pll - P12)/(R17 + RMF2) $ 

: W17' = (W17D - W17)/(TDELAY) $ 

: W18D = (PI 1 - P12)/(R18 + RVF2B) $ 

: W18' = (W18D - W18)/(TDELAY) $ 

: W19’ = 1/L19*(P10 - PLP - PL19) $ 

* pressure equations (psi) $ 

: PI' = 1/C1*(W1 - W7 - W8) $ 

: P2’ = 1/C2*(W7 - W4 - W6) $ 

: P3' = 1/C3*(W4 - W2 - W3) $ 

: P4' = 1/C4*(W2 + W3 - W5) $ 

: P5‘ = 1/C5*(W5 + W6 - W9) $ 

: P6' = 1/C6*(W8 - W10 - W15) $ 

: P7' = 1/C7*(W10 - Wll - W12) $ 

:P8'= 1/C8*(W11 +W12- W13)$ 
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P9' = 1/C9*(W9 + W13 - W14) $ 
P10’ = 1/C10*(W14 + W16 - W19) $ 
Pll'= 1/C11*(W15- W17-W18)$ 
P12' = 1/C12*(W17 + W18 - W16) $ 
PLP’ = 1/CACC*(W19 - Wl) $ 


* pressure sensor equation (transfer function of 1 ) (psi) $ 

: PHEAD = PI -P10$ 

* equation for the high pressure side of the pump (psi) $ 

: PHP = PLP + PUMP(W 1 ,PUMP\SPEED,PUMP\REF\SPEED) $ 


* pressure losses in the flow paths (psi) $ 

: PL1 = W1*R1 $ 

: PL4 = W4*R4 $ 

: PL5 = W5*R5 $ 

: PL7 = W7*R7 $ 

: PL8 = W8*R8 $ 

: PL9 = W9*R9 $ 

: PL10 = W10*R10 $ 

: PL11 = W11*R11 $ 

: PL12 = W12*R12 $ 

: PL14 = W14*R14 $ 

: PL15 = W15*R15 $ 

: PL16 = W16*R16 $ 

: PL19 = W19*R19 $ 

* include the data for the flow and pressure equations $ 
INCLUDE: "data.flow" $ 


* convert time in seconds to time in minutes for output purposesS 
: TIMEXMIN = TIME/60 $ 

ENDS 


**** defines **** 


* This file is located in /home/jacksme/marsyas/fum/defines/ on host edl4j. 

* It contains all the defined functions used in the SSFF TCS simulation $ 


* area of the pipe (ft A 2) $ 

DEFINE: AREA (DIA) = PI * (DIA/2)** * 2 $ 

* calculation of the reynolds number (unitless). $ 

DEFINE: RE (FLOW.VISC) = (FABS(FLOW) / (3600*AREA(DC))) * DC / (VISC * GC) $ 

* table for viscosity of fluid (lbf-sec/ft A 2) versus temperature (fah). 

obtained from the textbook entitled "fluid mechanics" by Streeter and 
wylie (p711). $ 

DEFINE: VISC (TEMP) = 

AF([ 0.0, 3.746E-5, 
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32.0, 3.746E-5, 

40.0, 3.229E-5, 

50.0, 2.735E-5, 

60.0, 2.359E-5, 

70.0, 2.050E-5, 

80.0, 1.799E-5, 

90.0, 1.595E-5, 

100.0, 1.424E-5, 

110.0, 1.284E-5, 

120.0, 1.168E-5, 

130.0, 1.069E-5, 

140.0, 0.98 IE-5, 

150.0, 0.905E-5 ], TEMP) $ 


* emperical formula for calculating the friction factor where the 

relative line roughness (k/dn) equals 0.00008 (number supplied by boeing). 
formula obtained from a paper entitled "a theoretical manual for g.l. von 
pragenau’s liquid seal simulation program" (may/june 1991) by dr. alan 
palazzolo. $ 

DEFINE: FRICTION (RE) = 0.0055*(l+(20000.0*0.00008+(10**6)/(RE+l))**(l/3)) $ 

* calculate kfactor given pressure drop @ a certain flow (unitless) $ 

DEFINE: KFACTOR (PRESSURE, FLOW) = PRESSURE/FABS(FLOW)**2 * 2 * RHOF 
* (AREA(DC))**2 * GC * (3600**2) * 144 $ 

* resistance (k-factor, flowrate) - calculate the resistance 
of a component (lbf-hr/lbm-in A 2). $ 

DEFINE: R (K, FLOW) = K*FABS(FLOW) / (2*RHOF*GC*144*(3600**2)*(AREA(DC))**2) $ 


* table defining valve angle versus loss factor supplied by allied-signal $ 
DEFINE: KVALVE (ANGLE) = 

AF([ 


0.0,30000000.0, 

2.8,: 

5995463.0, 

5.0, 

1 199092.0, 

6.0, 

599546.0, 

7.5, 

257805.0, 

10.0, 

101923.0, 

12.5, 

38971.0, 

15.0, 

21584.0, 

17.5, 

11691.0, 

20.0, 

8394.0, 

30.0, 

1559.0, 

40.0, 

539.6, 

50.0, 

251.0, 

60.0, 

125.0, 

75.0, 

47.0, 

90.0, 

22.0 ], ANGLE) $ 
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* table defining loss factor versus valve angle (degs) $ 
DEFINE: ANGLE (K) = 

AF([ 22.0, 90.0, 

47.0, 75.0, 

125.0, 60.0, 

251.0, 50.0, 

539.6, 40.0, 

1559.0, 30.0, 

8394.0, 20.0, 

11691.0. 17.5, 

21584.0, 15.0, 

38971.0, 12.5, 

101923.0, 10.0, 

257805.0, 7.5, 

599546.0, 6.0, 

1199092.0. 5.0, 

5995463.0, 2.8, 

30000000.0, 0.0 ], K) $ 


* this is a pump map (lbm/hr versus psi) obtained at "rspeed" rpm. 
rspeed/pspeed scales-up the flow to use a "rspeed" rmp map. 
(pspeed/rspeed)**2 scales-down the output pressure to a "pspeed" 
rmp map. $ 

DEFINE: PUMP (FLOW.PSPEED.RSPEED) = ((PSPEED/RSPEED)**2) * 
AF([ 

0.0,78.50, 

250.2. 78.50, 

500.5.78.50, 

1001.0, 78.50, 

1251.2,78.25, 

1501.5.78.00, 

1751.7.77.00, 

2002.0, 76.00, 

2252.2, 74.50, 

2502.5, 72.50, 

2752.7, 69.50, 

3003.0, 66.00, 

3253.2.61.50, 

3503.5, 54.00, 

3603.6, 50.00 ], FABS(FLOW)*(RSPEED/PSPEED)) $ 

* k-factor for a fluid line $ 

DEFINE: KLINE (FLOW,TEMP,LENGTH) = FRICTION((FABS(FLOW) / 
(3600*AREA(DI))) * DI / (VISC(TEMP) * GC))*LENGTH/DI $ 

* resistance for a fluid line $ 

DEFINE: RLINE (FLOW, TEMP, LENGTH) = 
KLINE(FABS(FLOW),TEMP,LENGTH)*FABS(FLOW) / 

(2*RHOF*GC* 144*(3600**2)*(AREA(DI))**2) $ 


21 



* convert watt to btu/sec $ 

DEFINE: QBTU(WATT) = WATT* * * * 3.4 12 1/3600 $ 


* temperature sensor time constant (hr) versus flow (lbm/hr) $ 
♦DEFINE: TC(FLOW) = 1/3600 * AF([ 29.96, 1 15.60, 
109.95, 68.03, 

230.44, 28.50, 

428.38, 12.94, 

812.40, 9.87, 

1224.14, 7.46], FLOW) $ 


* ua/ac vs mass flow for cp5 (btu/sec-f) $ 
♦DEFINE: UA56(FLOW) = 

AC56/3600*AF([ 0.0, 725.0, 

25.0, 725.0, 

50.0, 830.0, 

100.0, 915.0, 

300.0, 1070.0, 

500.0, 1070.0], FLOW) $ 


***** data-flow ***** 


* This file is located in /home/jacksme/marsyas/fum/data.flow on host edl4j. 

* It contains flow dynamic data and calculations for flow resistance 

* and inertance for the SSFF TCS simulation $ 

* note: all pressures are in units of lbf/in A 2 (psi) 

all flow rates are in units of lbm/hr $ 

* pump loss data $ 

EQUATION: PVPUMP = 16.7 $ 

: MDPUMP = 3000 $ 

: KPUMP = KFACTOR (P\PUMP, MDPUMP) $ 

* flow meter loss data $ 

: PFMETER = 1.1 $ 

: MDFMETER = 1300 $ 

: KFMETER = KFACTOR (PFMETER, MDFMETER) $ 

* shut-off valve loss data $ 

: PSOV = 0.2 $ 

: MDSOV = 3000 $ 

: KSOV = KFACTOR (PSOV, MDSOV) $ 

* temperature sensor loss data $ 

: PTEMP = 0.29 $ 

: MDTEMP = 600 $ 

: KTEMP = KFACTOR (PTEMP, MDTEMP) $ 

* quick disconnect loss data $ 

: PQD = 0.35 $ 

:MDQD= 1230$ 

: KQD = KFACTOR (PQD, MDQD) $ 

* pressure transducer loss data $ 
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: PPT = 1.29$ 

: MDPT = 600 $ 

: KPT = KF ACTOR (PPT, MDPT) $ 

* manual valve loss data $ 

* : PMAN = 0.25 $ 

* : MDMAN = 1230 $ 

* : KM AN = KF ACTOR (PMAN, MDMAN) $ 

* water-to-water (ww) heat exchanger loss data $ 

: PWW = 0.5 $ 

: MDWW = 1100$ 

: KWW = KF ACTOR (PWW, MDWW) $ 

* water-to-air (wa) heat exchanger loss data $ 

: PWA = 1.0$ 

: MDWA = 180 $ 

: KWA = KF ACTOR (PWA, MDWA) $ 

* coldplate type cp- 1-2-3 loss data $ 

: PCP123 = 0.6 $ 

:MDCP 123 = 1000.0$ 

: KCP123 = KFACTOR (PCP123, MDCP123) $ 

* coldplate type cp-4-5 -6 loss data $ 

: PCP456 = 0.6 $ 

: MDCP456 = 300.0 $ 

: KCP456 = KFACTOR (PCP456, MDCP456) $ 

* coldplate type cp-7-8 loss data $ 

: PCP78 = 1.0 $ 

: MDCP78 = 300.0 $ 

: KCP78 = KFACTOR (PCP78, MDCP78) $ 

* furnace jacket loss data $ 

: PFURN = 0.7 $ 

: MDFURN = 150 $ 

: KFURN = KFACTOR (PFURN, MDFURN) $ 

* 3 foot flex line (1/2" dia) loss data $ 

: PFLEX = 0.77 $ 

: MDFLEX = 1230 $ 

: KFLEX = KFACTOR (PFLEX, MDFLEX) $ 

* 90 degree pipe bend loss data $ 

* : P90B = 0.038 $ 

* : MD90B = 3000 $ 

* : K90B = KFACTOR (P90B, MD90B) $ 

* resistance calculation for the flow control valves $ 

: RVF1A = R( K V AL VE(EM 1 AW AL VE\ANGLE) , W 2) $ 

: RVC = R(KVALVE(CORE\VALVE\ANGLE),W 1 3) $ 

: RVF2B = R(KV ALVE(EM2B\V ALVE\ANGLE) ,W 1 8) $ 

* resistance calculation for the manual valves (static paths) $ 

: RMF1 = R(KV ALVE( AMF 1 ),W3) $ 

: RMC = R(KVALVE(AMC),W6) $ 

: RMF2 = R(KVALVE(AMF2),W17) $ 

* lumped pipe line lengths (ft) for each path $ 

: LL1 = 12.0/12 $ 
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: LL2 = 120.0/12$ 

:LL3 = 72.0/12 $ 

: LL4 = 50.0/12 $ 

:LL5 = 52.0/12$ 

: LL6 = 72.0/12 $ 

: LL7 = 20.0/12 $ 

: LL8 = 20.0/12 $ 

:LL9 = 40.0/12$ 

:LL10 = 24.0/12$ 

:LL11 =24.0/12$ 

: LL12 = 24.0/12 $ 

: LL13 = 24.0/12 $ 

:LL14 = 20.0/12$ 

:LL15 = 26.5/12$ 

:LL16 = 42.0/12$ 

: LL17 = 72.0/12$ 

:LL18 = 72.5/12$ 

: LL19 = 146.0/12$ 

* lumped resistance calculations for each flow path (lbf-hr/lbm-in A 2) $ 

: R1 = R(KPUMP,Wl)/2 + R(KFLEX,W1) + R(KQD,W1) + R(KPT,W1) 
+RLINE(W1, T, LL1) $ 

: R2 = R(KFURN,W2) + 6*R(KQD,W2) + R(KFMETER,W2) + R(KTEMP,W2) 
+R(KSOV,W2) 

+ 2*R(KFLEX,W2) + RLINE(W2, T, LL2) $ 

: R3 = R(KCP78,W3) + R(KCP456,W3) + 2*R(KCP78,W3) + R(KSOV,W3) 

+ 12*R(KQD,W3) + 4*R(KFLEX,W3) + RLINE(W3, T, LL3) $ 

: R4 = RLINE(W4, T, LL4) $ 

: R5 = RLINE(W5, T, LL5) $ 

: R6 = 4*R(KCP456,W6) + R(KCP123,W6) + R(KCP78,W6) + 14*R(KQD,W6) 

+ 2*R(KFLEX,W6) + RLINE(W6, T, LL6) $ 

: R7 = RLINE(W7, T. LL7) $ 

: R8 = RLINE(W8, T, LL8) $ 

: R9 = RLINE(W9, T, LL9) $ 

: RIO = R(KWA,W10) + 2*R(KQD,W10) + RLINE(W10, T, LL10) $ 

: R 1 1 = R(KCP123,W1 1) + RLINE(W1 1, T, LL1 1) $ 

: R12 = R(KCP123,W12) + RLINE(W12, T, LL12) $ 

: R13 = R(KCP78,W13) + 2*R(KFLEX,W13) + R(KTEMP,W13)+4*R(KQD,W13) 
+ R(KFMETER,W13) + RLINE(W13, T, LL13) $ 

: R14 = RLINE(W14, T, LL14) $ 

: R15 = RLINE(W15, T, LL15) $ 

: R16 = RLINE(W16, T, LL16) $ 

: R17 = R(KCP78,W17) + R(KCP456,W17) + 2*R(KCP78,W17) + R(KSOV,W17) 
+ 12*R(KQD,W17) + 4*R(KFLEX,W17) + RLINE(W17, T, LL17) $ 

: R18 = R(KFURN,W18) + 6*R(KQD,W18) + R(KFMETER,W 1 8) 
+R(KTEMP,W18) 

+ R(KSOV,W18) + 2*R(KFLEX,W 1 8) + RLINE(W18, T, LL18) $ 

: R19 = R(KPUMP,W19)/2 + R(KWW,W19) + R(KPUMP,W19)/2 + R(KQD,W19) 
+ 3*R(KFLEX,W 19) + R(KPT,W19) + RL1NE(W19, T, LL19) $ 
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* line cross sectional area (ft A 2) $ 

: A = AREA(DI) $ 

* flex line length (ft) $ 

: FLEX = 3.0 $ 

* lumped line inertances (hr-sec-lbf/lbm-in A 2) plus equivalent rack inertances 
calculations for each flow path $ 

: LI = (LLl+FLEX)/(3600* 144*GC*A) $ 

: L4 = LL4/(3600*144*GC*A) + LEQ1/2 $ 

: L5 = LL5/(3600* 144*GC*A) + LEQ1/2 $ 

: L7 = LL7/(3600* 144*GC*A) + L6/2 $ 

: L8 = LL8/(3600* 144*GC*A) $ 

: L9 = LL9/(3600* 144*GC*A) + L6/2 $ 

: L10 = LL 1 0/(3600* 144*GC* A) $ 

: Lll = LL1 1/(3600* 144*GC*A) + L13/2 $ 

: L12 = LL12/(3600*144*GC*A) + L13/2 $ 

: L14 = LL14/(3600* 144*GC*A) $ 

: L15 = LL15/(3600*144*GC*A) + LEQ2/2 $ 

: L16 = LL16/(3600* 144*GC*A) + LEQ2/2 $ 

: L19 = (LL19+2*FLEX)/(3600*144*GC*A) $ 

: L2 = (LL2+2*FLEX)/(3600* 1 44*GC * A) $ 

: L3 = (LL3+4*FLEX)/(3600*144*GC*A) $ 

: L6 = (LL6+2*FLEX)/(3600*144*GC*A) $ 

: L13 = (LL13+2*FLEX)/(3600*144*GC*A) $ 

: L17 = (LL17+2*FLEX)/(3600*144*GC*A) $ 

: L18 = (LL 1 8+2*FLEX)/(3600* 1 44*GC* A) $ 

: LEQ1 = L2*L3/(L2+L3) $ 

: LEQ2 = L17*L18/(L17+L18) $ 

* node capacitances (lbf-hr/lbm-in A 2-sec) $ 

: Cl =0.1 $ 

: C2 = 0. 1 $ 

: C3 = 0.1 $ 

: C4 = 0.1 $ 

: C5 = 0.1 $ 

: C6 = 0.1 $ 

: C7 = 0.1 $ 

: C8 = 0.1 $ 

: C9 = 0.1 $ 

: CIO = 0.1 $ 

: C 1 1 =0.1 $ 

: C12 = 0.1 $ 

: CACC = 0.1 $ 


**** £ 

* This file is located in /home/jacksme/marsyas/furn/f_act on host edl4j. 

* It contains the actuator/valve models for each flow control path $ 


MODEL: EM 1 AW AL VE\ACTU ATOR $ 



INPUTS: VA $ 

OUTPUTS: AV $ 

* wv - valve speed response to the input speed voltage command (va) $ 

* tau - delay between commanded and actual valve speed response (secs) $ 

* kv - voltage to speed conversion constant (degs/sec/volts) $ 

* av - valve position (degs) $ 

EQUATION: WV' = 1/TAU* * * (KV*VA - WV) $ 

: AV’ = WV $ 

: KV = 1.0 $ 

: TAU = 0.01 $ 

END $ 

MODEL: COREWALVEVACTUATOR $ 

INPUTS: VA $ 

OUTPUTS: AV $ 

* wv - valve speed response to the input speed voltage command (va) $ 

* tau - delay between commanded and actual valve speed response (secs) $ 

* kv - voltage to speed conversion constant (degs/sec/volts) $ 

* av - valve position (degs) $ 

EQUATION: WV' = 1/TAU*(KV*VA - WV) $ 

: AV' = WV $ 

: KV= 1.0$ 

: TAU = 0.01 $ 

ENDS 

MODEL: EM2BW ALVEYACTU ATOR $ 

INPUTS: VA $ 

OUTPUTS: AV $ 

* wv - valve speed response to the input speed voltage command (va) $ 

* tau - delay between commanded and actual valve speed response (secs) $ 

* kv - voltage to speed conversion constant (degs/sec/volts) $ 

* av - valve position (degs) $ 

EQUATION: WV' = 1/TAU*(KV*VA - WV) $ 

: AV' = WV $ 

:KV= 1.0$ 

: TAU = 0.01 $ 

END $ 

**** f cont **** 


* This file is located in /home/jacksme/marsyas/fum/f_cont on host edl4j. 

* This file contains the submodels for the continuous and discrete 

* flow control algorithms. First, the controllers were designed in the 

* discrete domain and then converted to the continuous domain using 

* pole-zero mapping (ref. Discrete Time Control Systems by Ogata). $ 

* emla discrete flow control submodel $ 

MODEL: EMlA\FLOW\CONTROLLER\DISCRETE,DISCRETE $ 

INPUTS: W2.W2SET $ 

OUTPUTS: EMI A\VALVE\CMD $ 
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* bumpless transfer function - digital low pass filter that does not 
allow the set point to change rapidly $ 

EQUATION: BSET = 0.90*BSETOLD + 0.10*W2SET $ 

: BSETOLD% = BSET $ 

* flow error based on filtered set-point $ 

: ERR = BSET - W2 $ 

* digital low pass noise filter does not allow the current error 
to effect the computed output error by more than 10 percent, 
the previously computed output error comprises 90 percent of 
the computed output error $ 

: FERR = 0.90*FERROLD + 0. 10*ERR $ 

: FERROLD% = FERR $ 

* compute error derivative $ 

: FEDOT = (FERR - FERROLD)/1.0 $ 

* velocity pd control algorithm $ 

: U = FERR*0.015 + FEDOT*0.1 $ 

* valve command limiter $ 

: VCMD = 5.0 IF U.GT.5.0 ELSE 
-5.0 IF U.LT.-5.0 ELSE 
U$ 

: EM 1 AW ALVE\CMD = VCMD $ 

ENDS 

* emla continuous flow control submodel $ 

MODEL: EMlA\FLOW\CONTROLLER\CONTINUOUS $ 

INPUTS: W2.W2SET $ 

OUTPUTS: EM 1 AW ALVE\CMD $ 

* bumpless transfer function - digital low pass filter that does not 
allow the set point to change rapidly $ 

EQUATION: BW2SET' = -.1054*BW2SET + .1054*W2SET $ 

* flow error based on filtered set-point $ 

: ERR = BW2SET - W2 $ 

* digital low pass noise filter does not allow the current error 
to effect the computed output error by more than 10 percent, 
the previously computed output error comprises 90 percent of 
the computed output error $ 

: FERR = TF( 1 ,0,0. 1 054, 1 .0,0. 1054.ERR) $ 

* velocity pd control algorithm $ 

: U1A = TF( 1 1 ,.0 1 5,TDELAY, 1 .0,FERR) $ 

* valve command limiter $ 

: VCMD 1 A = 5.0 IF U1A.GT.5.0 ELSE 
-5.0 IF U1A.LT.-5.0 ELSE 
U1A $ 

* computer delay (= 1.0 sec) $ 

: EM 1 AW ALVE\CMD' = -2*EM1A\VALVE\CMD + 2*VCMD1A $ 
ENDS 

* core discrete flow control submodel $ 

MODEL: CORE\FLOW\CONTROLLER\DISCRETE,DISCRETE $ 
INPUTS: W13.W13SETS 
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OUTPUTS: CORE\VALVE\CMD $ 

* bumpless transfer function - digital low pass filter that does not 
allow the set point to change rapidly $ 

EQUATION: BSET = 0.90*BSETOLD + 0.10*W13SET $ 

: BSETOLD% = BSET $ 

* flow error based on filtered set-point $ 

: ERR = BSET -W13$ 

* digital low pass noise filter does not allow the current error 
to effect the computed output error by more than 10 percent, 
the previously computed output error comprises 90 percent of 
the computed output error $ 

: FERR = 0.90*FERROLD + 0.10*ERR $ 

: FERROLD% = FERR $ 

* compute error derivative $ 

: FEDOT = (FERR - FERROLD)/1.0 $ 

* velocity pd control algorithm $ 

: U = FERR*0.015 + FEDOT*0.1 $ 

* valve command limiter $ 

: VCMD = 5.0 IF U.GT.5.0 ELSE 
-5.0 IF U.LT.-5.0 ELSE 
U$ 

: CORE\VALVE\CMD = VCMD $ 

ENDS 

* core continuous flow control submodel $ 

MODEL: CORE\FLOW\CONTROLLER\CONTINUOUS $ 

INPUTS: W13,W BSET $ 

OUTPUTS: CORE\VALVE\CMD $ 

* bumpless transfer function - digital low pass filter that does not 
allow the set point to change rapidly $ 

EQUATION: BW13SET' = -,1054*BW13SET + ,1054*W1 3SET $ 

* flow error based on filtered set-point $ 

: ERR = BW13SET - W13 $ 

* digital low pass noise filter does not allow the current error 
to effect the computed output error by more than 10 percent, 
the previously computed output error comprises 90 percent of 
the computed output error $ 

: FERR = TF( 1,0,0. 1054, 1.0,0. 1054, ERR) $ 

* velocity pd control algorithm $ 

: U = TF( 1 ,. 1 ,.0 1 5,TDEL A Y, 1 .0.FERR) $ 

* valve command limiter $ 

: VCMD = 5.0 IF U.GT.5.0 ELSE 
-5.0 IF U.LT.-5.0 ELSE 
US 

* computer delay (= 1.0 sec) $ 

: CORE\VALVE\CMD' = -2*CORE\VALVE\CMD + 2*VCMD $ 
ENDS 

* em2b discrete flow control submodel $ 

MODEL: EM2B\FLOW\CONTROLLER\DISCRETE,DISCRETE $ 


28 



INPUTS: W18.W18SETS 
OUTPUTS: EM2BYV ALVE\CMD $ 

* bumpless transfer function - digital low pass filter that does not 
allow the set point to change rapidly $ 

EQUATION: BSET = 0.90*BSETOLD + 0.10*W18SET $ 

: BSETOLD% = BSET $ 

* flow error based on filtered set-point $ 

: ERR = BSET - W18$ 

* digital low pass noise filter does not allow the current error 
to effect the computed output error by more than 10 percent, 
the previously computed output error comprises 90 percent of 
the computed output error $ 

: FERR = 0.90*FERROLD + 0.10*ERR $ 

: FERROLD% = FERR $ 

* compute error derivative $ 

: FEDOT = (FERR - FERROLD) / 1.0 $ 

* velocity pd control algorithm $ 

: U = FERR*0.015 + FEDOT*0.1 $ 

* valve command limiter $ 

: VCMD = 5.0 IF U.GT.5.0 ELSE 
-5.0 IF U.LT.-5.0 ELSE 
U $ 

: EM2BW ALV E\C M D = VCMD $ 

END $ 

* em2b continuous flow control submodel $ 

MODEL: EM2B\FLOW\CONTROLLER\CONTINUOUS $ 

INPUTS: W18,W18SET$ 

OUTPUTS: EM2B\VALVE\CMD $ 

* bumpless transfer function - digital low pass filter that does not 
allow the set point to change rapidly $ 

EQUATION: BW18SET' = -.1054*BW18SET + .1054*W18SET $ 

* flow error based on filtered set-point $ 

: ERR = BW18SET - W18 $ 

* digital low pass noise filter does not allow the current error 
to effect the computed output error by more than 10 percent, 
the previously computed output error comprises 90 percent of 
the computed output error $ 

: FERR = TF(1,0, 0.1054, 1.0, 0.1054, ERR) $ 

* velocity pd control algorithm $ 

: U = TF( 1,1 ,.0 1 5,TDELAY, 1 .0,FERR) $ 

* valve command limiter $ 

: VCMD = 5.0 IF U.GT.5.0 ELSE 
-5.0 IF U.LT.-5.0 ELSE 
U$ 

* computer delay (= 1.0 sec) $ 

: EM2BYV ALVE\CMD' = -2*EM2B\VALVE\CMD + 2*VCMD $ 
ENDS 

**** p_cont **** 
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* This file is located in /horne/jacksme/marsyas/furn/p_cont on host edl4j. 

* It contains the continuous and discrete pressure control algorithms $ 

* continuous pressure controller $ 

MODEL: PRESSURE\CONTROLLER\CONTINUOUS $ 

INPUTS: PHEAD.PHEADSET $ 

OUTPUTS: PUMP\SPEED\CMD $ 

EQUATION : PERR = PHEADSET - PHEAD $ 

* controller transfer function = 

* ko*(s A 2 + (br2+br3)*s + br2*br3) / (s A 2 + (brl+br4)*s + brl*br4) $ 
EQUATION : XI* = -A1*X1 + X2 + (KO*Bl - Al*KO)*PERR $ 

: X2' = KO*BO*PERR $ 

: CMD = (XI + KO*PERR)*2*PI/60 $ 

* controller coefficient definitions $ 

: A1 = BR1+BR4 $ 

: B1 = BR2+BR3 $ 

: BO = BR2*BR3 $ 

* controller gain (volts/psi) $ 

: KO = 6000 $ 

* controller break frequencies (rads/sec) $ 

:BR1 =0.0$ 

: BR2 = 0.1*2*PI $ 

: BR3 = 1.0*2*PI $ 

: BR4 = 10.0*2*PI $ 

* computer delay (= 1.0 sec) $ 

: PUMP\SPEED\CMD' = -2*PUMP\SPEED\CMD + 2*CMD $ 

ENDS 

* discrete pressure controller $ 

MODEL: PRESSUREVCONTROLLERYDISCRETE, DISCRETE $ 

INPUTS: PHEAD.PHEADSET $ 

OUTPUTS: PUMP\SPEED\CMD $ 

EQUATION: E = PHEADSET - PHEAD $ 

* difference equation obtained by converting the continuous pressure 
controller to the discrete domain using pole-zero mapping (ref. Discrete 
Time Control Systems by Ogata) $ 

: U = UD + KSTAR*(E - (EXP(-BR2*T)+EXP(-BR3*T))*ED% 

+ EXP(-BR2*T)*EXP(-BR3*T)*ED) $ 

: ED%% = E $ 

: UD% = U $ 

* discrete controller gain (volts/psi) $ 

: KSTAR = KO*(BR2+BR3)/(BR4*(2-EXP(-BR2*T)-EXP(-BR3*T))) $ 

* continuous controller gain (volts/psi) $ 

: KO = 6000 $ 

* continuous controller break frequencies (rads/sec) $ 

: BR2 = 0.1*2*PI $ 

: BR3 = 1.0*2*PI $ 

: BR4 = 10.0*2*PI $ 

:T = 1.0$ 
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* convert pump speed from rev/min to rad/sec $ 

: PUMP\SPEED\CMD = U*2*PI/60 $ 

END $ 

**** pump_sc **** 

* This file is located in /home/jacksme/marsyas/furn/pump_sc on host edl4j. 

* It contains the pump speed controller model supplied by Allied Signal 

* (pump vendor). $ 

MODEL: PUMP\SPEED\CONTROLLER $ 

INPUTS: PUMP\SPEED\CMD,PUMP\SPEED\SEN $ 

OUTPUTS: PUMPWOL'RCMD $ 

* speed command limiter in (radians/second) $ 

EQUATION: SPEED\CMD\LIM = 22222.0*(2*PI/60) IF PUMP\SPEED\CMD.GT.22222.0*(2*PI/60) 
ELSE 

0.0 IF PUMP\SPEED\CMD.LT.0.0 ELSE 
PUMP\SPEED\CMD $ 

* ^pnQPfl <inppd filtpr $ 

: SPEEDVFBK’ = (PUMP\SPEED\SEN - SPEED\FBK)/0.00796 $ 

* calculate speed error (radians/second) $ 

: SPEEDXERR = (SPEED\CMD\LIM - SPEED\FBK)/2327.1 $ 

* speed controller $ 

* : VOLT = 120*TF( 1 ,. 1 24, 1 ,.3,0,SPEED\ERR) $ 

: XT = SPEED\ERR $ 

: VOLT = 120.0/.3*(.124*SPEED\ERR + XI) $ 

* voltage limiter $ 

: PUMPYV OLT\CMD = 120 IF VOLT.GT.120 ELSE 
0 IF VOLT.LT.O ELSE 
VOLTS 

ENDS 

**** pump_motor **** 

* This file is located in /home/jacksme/marsyas/fum/pump_motor on host edl4j. 

* It contains the pump motor and speed sensor model supplied by Applied 

* Signal (pump vendor). $ 

MODEL: PUMPVMOTOR $ 

INPUTS: PUMPW OLTVCMD $ 

OUTPUTS: PUMP\SPEED,PUMP\SPEED\SEN,PUMP\CURR\SEN $ 

EQUATION: VOLTVERR = PUMP\VOLT\CMD - VEMF $ 

* back emf gain (volt/rads/sec) $ 

:KE = 0.04185$ 

* back emf calculation (volts) $ 

: VEMF = KE*PUMP\SPEED\SEN $ 

* lm - motor inductance (henries) $ 

:LM = 0.0012$ 

* rm - motor resistance (ohms) $ 

:RM = 0.618$ 
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* motor impedance transfer function $ 

* : PUMP\CURR\SEN = TF(1,0.0,1/RM,LM/RM,1.0,VOLT\ERR) $ 

: PUMP\CURR\SEN' = -RM/LM*PUMP\CURR\SEN + l/LM*VOLTVERR $ 

* torque gain (in-lb/amp) $ 

: KT = 0.33125$ 

* motor torque calculation (in-lb) $ 

: TM = KT*PUMP\CURR\SEN $ 

* actual torque calculation (in-lb) $ 

: TA = TM - TL $ 

* empirical torque load calculation (in-lb) $ 

: TL = 1.9875*(PUMP\SPEED\SEN)**2/1979**2 $ 

* rotor and load inertia (in-lb-s A 2) $ 

: JML = 500E-6/16 $ 

* derivative of angular speed equals torque/inertia (rads/sec) $ 

: PUMP\SPEED\SEN' = TA/JML $ 

* output equation to pump map $ 

: PUMP\SPEED = PUMP\SPEED\SEN*60/(2*PI) $ 

END $ 

**** **** 


* This file is located in /home/jacksme/marsyas/furn/sim on host edl4j. 

* This file contains the simulation module for the main SSFF TCS model. 

* It includes the input excite functions, the initialization statements, 

* the discrete submodel schedules, etc. $ 

SIMULATE: MAIN $ 

* excitation of the inputs in the main model $ 

EXCITE: AMF1 =2.13071 1E+01 $ 

: AMC = 2.1 2393 1E+01 $ 

: AMF2 = 2.123383E+01 $ 

: PUMP\REF\SPEED = 18900.0 $ 

: T = 80.0 $ 

: PHEADSET = 10.0 $ 

: W2SET = 260.0 $ 

: W13SET = 386.0$ 

:W18SET = 260.0$ 

* set initial values for the states in the main model $ 

INITIALIZE: MAIN, PLP = 2.600000E+01, 

W1 = 1.056000E+03, 

W2 = 2.600000E+02, 

W3 = 5.000000E+01, 

W4 = 3. 100000E+02, 

W5 = 3.100000E+02, 

W6 = 5.000000E+01, 

W7 = 3.600000E+02, 

W8 = 6.960000E+02, 

W9 = 3.600000E+02, 

W10 = 3.860000E+02, 
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Wll = 193.0, 

W12 = 193.0, 

W13 = 3.860000E+02, 

W14 = 7. 460000 E+02, 

W15 = 3.100000E+02, 

W16 = 3.100000E+02, 

W17 = 5.000000E+01, 

W18 = 2.600000E+02, 

W19 = 1.056000E+03, 

PI = 4.513639E+01, 

P2 = 4.512252E+01, 

P3 = 4.509572E+0 1 , 

P4 = 3.524069E+01, 

P5 = 3.521281E+01, 

P6 = 4.5093 19E+01, 

P7 = 4.040686E+01, 

P8 = 4.037880E+01, 

P9 = 3.518509E+01, 

P10 = 3.513639E+01, 

PI 1 = 4.507899E+01, 

P12 = 3.515890E+01 $ 

* set initial values in the pump speed controller submodel $ 

INITIALIZE: PUMP\SPEED\CONTROLLER, SPEEDVFBK = 18900.0*5.717402E-01*(2*PI/60) $ 
: PUMP\SPEED\CONTROLLER, XI = 1.214235E-01 $ 

* set initial values in the pump motor submodel $ 

INITIALIZE: PUMP\MOTOR, PUMP\CURR\SEN = 1.961725E+00 $ 

: PUMPVMOTOR, PUMP\SPEED\SEN = 1 8900.0*5.7 17402E-01*(2*PI/60) $ 

* set initial values in the valve actuator submodel $ 

INITIALIZE: EM 1 AW ALVEVACTU ATOR, WV = 0.0, AV = 5.310246E+01 $ 

: COREW AL VEYACTU ATOR, WV = 0.0, AV = 7.984484E+01 $ 

: EM2B\ V AL V E\ ACTU ATOR , WV = 0.0, AV = 5.291287E+01 $ 

* set initial values in the continuous pressure controller submodel $ 

INITIALIZE: PRESSURE\CONTROLLER\CONTINUOUS, 

Xl=18900.0*5.7 17402E-01, 

X2= 1 8900.0*5.7 17402E-0 1 * 10*2*PI, 

PUMP\SPEED\CMD = (18900.0*5.717402E-01)*2*PI/60 $ 

* set initial values in the discrete pressure controller submodel $ 

INITIALIZE: PRESSURE\CONTROLLER\DISCRETE, UD = 18900.0*5.717402E-01 $ 

* set initial value in the continuous and discrete flow controllers $ 

INITIALIZE: EMlA\FLOW\CONTROLLER\CONTINUOUS, BW2SET = 260.0 $ 

INITIALIZE: CORE\FLOW\CONTROLLER\CONTINUOUS, BW13SET = 386.0 $ 

INITIALIZE: EM2B\FLOW\CONTROLLER\CONTINUOUS, BW18SET = 260.0 $ 

INITIALIZE: EM 1 A\FLOW\CONTROLLER\DISCRETE, BSETOLD = 260.0 $ 
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INITIALIZE: COREVFLOVACONTROLLERYDISCRETE, BSETOLD = 386.0 $ 

INITIALIZE: EM2B\FLOW\CONTROLLER\DISCRETE, BSETOLD = 260.0 $ 

* numerically integrate using runga-kutta 4th order method $ 

INTEGRATE: RK4, TIMESTEP, TSTEP $ 

* set simulation stop time $ 

STOPIF : TIME.GE.STOPT $ 

* plot the main outputs every 100 integration steps $ 

PLOT : LINE AR(0,STOPT, 1 00),\OUTPUTS\ $ 

* schedule each discrete controller every 1.0 second $ 

SCHEDULE: PRESSURE\CONTROLLER\DISCRETE,FMOD(TIME+TSTEP/2, 1 .0).LT TSTEP $ 
SCHEDULE: EM 1 A\FLOW\CONTROLLER\DISCRETE,FMOD(TIME+TSTEP/2, 1 .0).LT TSTEP $ 
SCHEDULE: CORE\FLOW\CONTROLLER\DISCRETE,FMOD(TIME+TSTEP/2, 1 .0).LT.TSTEP $ 
SCHEDULE: EM2B\FLOW\CONTROLLER\DISCRETE,FMOD(TIME+TSTEP/2, 1 .0).LT. TSTEP $ 

END $ 
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